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We study the parametric amplification of super-Hubble-scale scalar metric fluctuations at the end 
of inflation in some specific two-field models of inflation, a class of which is motivated by hybrid 
inflation. We demonstrate that there can indeed be a large growth of fluctuations due to parametric 
resonance and that this effect is not taken into account by the conventional theory of isocurvature 
perturbations. Scalar field interactions play a crucial role in this analysis. We discuss the conditions 
under which there can be nontrivial parametric resonance effects on large scales. 



PACS numbers: 98.80Cq 

I. INTRODUCTION 

It was recently suggested that parametric resonance 
during the reheating phase of an inflationary Universe 
may lead to an exponential amplification of super- 
Hubble scale gravitational fluctuations. If true, this 
would affect the usual predictions of inflationary mod- 
els for observables such as the matter power spectrum 
and the spectrum of cosmic microwave anisotropies. In 
particular, it would require the coupling constants in the 
particle physics model of inflation to be exponentially 
smaller than previously thought in order that the theory 
does not generate a too large amplitude for the fluctua- 
tions. 

In Ref. || it was shown that, although there are no 
causality constraints which prohibit the amplification of 
super-Hubble (but sub-horizon) modes during reheat- 
ing, the effect does not occur in a simple massive scalar 
field model of chaotic inflation based on the potential 
V{4>) = m 2 cj} 2 /2 (Here, is the inflaton field). This is 
true even beyond the linear analysis H . Similarly, there 
is no effect for a quartic potential 0] V{4>) = A<ji> 4 /4, 
nor for a 
tic terms 
in Refs 



otential containing both quadratic and quar- 
I . These results agree with the earlier analyses 
^ and 0. 

It was then suggested jll]] that the amplification of 
super-Hubble-scale modes would occur for two field mod- 
els of inflation, e.g. for a model with potential 



V(<t>, X ) = ^mV + ^W- 



(1) 



where as before <f> is the inflaton field and x is a second 
scalar matter field. This model had earlier been analyzed 
by Taruya and Nambu who claimed that the isocur- 
vature mode of the fluctuations will be parametrically 
amplified during reheating. However, as was shown in 
Refs. Jll| and (and more recently in Jl5|]), the fluctu- 
ations in the \ field are exponentially suppressed during 
inflation for values of the coupling constant for which the 



equation of motion of the metric fluctuations corresponds 
to broad resonance, thus rendering the effect studied in 
Ref. [|ll| completely inefficient. 

The suppression of fluctuations in the \ field which 
renders the parametric amplification of gravitational fluc- 
tuations ineffective in the model given by (j|) occurs since 
during inflation the induced mass m x of the x field which 
is given by m x — g\<j>\ is larger than the Hubble expan- 
sion parameter H, and hence, as can be easily seen by 
considering the equation of motion 



,k 



8 x + W5 x +{- + g 2 4> 2 )8x = 



(2) 



for the linearized fluctuation of the x field with comoving 
wave number k (the scale factor is denoted by a(t)), 5\ 
undergoes damped oscillatory motion. 

A model in which m x < H during the stage of inflation 
when scales of cosmological interest today exit the Hub- 
ble radius was recently studied by Bassett and Viniegra 
p6[ . It is a two field model given by the potential 



Vfax) = \^ + \g 2 4> 2 x 2 - 



(3) 



In the absence of metric fluctuations, this model was 
studied in detail in Ref. (see also jl^|), where it was 
shown that for values of the coupling constants satisfying 



g 



2X 



(4) 



long wavelength modes (k ~ 0) are in the first broad 
instability band of the Floquet-type equation of motion 
derived from (||) after field rescaling which describes the 
parametric resonance of matter fluctuations in an unper- 
turbed expanding space-time. Bassett and Viniegra |l~6| ] 
showed that in this model the quantity £ increases 
exponentially during the initial stages of reheating. Note 
that C is a measure of the curvature fluctuations and is 
believed to be conserved on super-Hubble scales in the 
absence of isocurvature fluctuations (see e.g. pOl for a 
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review of the theory of cosmological fluctuations). How- 
ever, since the model given by @ admits isocurvature 
fluctuations, a growth of £ on super-Hubble modes is ex- 
pected also in the "usual" analysis of the evolution of 
fluctuations in inflationary cosmology. 

In this paper we take a closer look at the theory given 
by the potential (||). Subject to certain assumptions on 
initial conditions of the background field dynamics we re- 
cover results similar to Bassett and Viniegra: exponential 
growth of £ during the initial stages of reheating. Fur- 
thermore, we demonstrate that this effect is not taken 
into account by the conventional theory of isocurvature 
perturbations. We then discuss some criteria which an 
inflationary Universe model must satisfy in order to have 
substantial parametric growth of £ during reheating via 
scalar field interactions. We argue that these conditions 
are naturally satisfied in some models of hybrid inflation, 
and we study a couple of concrete examples in which 
super-Hubble-scale gravitational fluctuations grow expo- 
nentially during reheating (some other examples where 
exponential growth of super-Hubble-scale modes could 
occur are given in pf[). 



II. MASSLESS TWO FIELD MODEL 
RECONSIDERED 

In this section we will consider the two-field model with 
potential (^) in which Bassett and Viniegra jl6| showed 
the parametric resonance of super-Hubble scale gravita- 
tional modes. This model has been studied in detail in 
[ 17 1 in the absence of gravitational perturbations (see also 

[!)•. 

As in our previous paper [pj, we shall work in longi- 
tudinal gauge in which the metric including linearized 
scalar metric fluctuations takes on the form 



ds 2 



(I - 2$)dt 2 - a 2 {t)(l + 2V)dx i dx i 



(5) 



where the x l are the comoving spatial coordinates and t is 
physical time. Since for the matter model considered the 
off-diagonal components of the spatial part of the energy- 
momentum tensor vanish, the corresponding components 
of the Einstein equations imply |^0| \& = For the 
moment we will write down the equations for a general 
system with multiple scalar fields (i = 1, n), and only 
at a later stage will we specialize to the specific model 
considered in this section. 

The remaining independent equations of motion for lin- 
earized perturbations in this Einstein-Higgs system are 
the perturbed energy constraint and momentum con- 
straint equations as well as the equations of motion for 
the Higgs field perturbations 5 fa: 



3#$ 



m 2 $ 



a~ 

AttG / I.-,.., 

i=l 



(6) 



n 

■n 



$ + H§ = 4ttG <i>i 



(7) 



ZHSfa 



k 2 
7^' 



- 2V n $ , 



(8) 



where k denotes the comoving wavenumber, V,i indicates 
the derivative of V with respect to fa, and G is Newton's 
constant. 

The Sasaki-Mukhanov J2^] variables for the n matter 
fields are 



Qi = Sfa + ^$ (9) 
ti 

and satisfy the following system of equations: 



i,2 " o f 3 
Qi + 3HQi + —Qt + y\V tij Hs<t>ih)-\Qi = ■ 



(10) 



We will now specialize to our two field model. If the 
homogeneous part of the second scalar field x vanishes 
(which will not be true if parametric resonance is to excite 
gravitational fluctuations - see later) , then the inflaton <f> 
during the initial stages of reheating (when back-reaction 
effects are negligible) oscillates as follows [|l7j 



far]) = a l facn(x - x , , 



(11) 



where fa is the amplitude of (f> at the end of the slow- 
rolling period, an is the Jacobi elliptic cosine function, 
r\ is conformal time and x = V\fa is a rescaled dimen- 
sionless conformal time coordinate. Following |l7| , it is 
insightful to rescale all the fields / by the scale factor 
and use new fields f = af. The equations of motion ( |Io| ) 
then become 



[k 2 (- 3cn 2 (x — xq, —7=)] 1 

V2 



a 

g 2 1 
-2^— cn(x - x , —=)xQx + M<t>4>Q<p + M^C 



(12) 



A0 O 
2 a , 9 



V2' 



1 



Qx + W + \cn 2 {x-x a ,-=)]Q 



.'/" 1 



(13) 



= ~ 2 l^~ cn ( x ~ x °' ^ )xQ<l> + M 4>xQ4> + M xxQx 

where, following the notation of (TgJ (up to a prefactor 
containing a), we have have used the abbreviation 

M Ma = — [-ff^i --Mi<h ~-<h)] > (14) 

Qj Q/J.J. Qj Qi 

and where 
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\4>l 



(15) 



H V 



In the model considered, the time-averaged equation of 
state is that of radiation and hence the scale factor is 
linear in rj and thus a =0. 

Neglecting for a moment the terms on the right hand 
side of the equations, both equations ( fl2|) and (|l3| ) are of 
the form of Lame equations. Lame equations in the con- 
text of reheating were first noticed in || and then studied 
in detail in fEj , (l8| and The coefficient in front of 
the cn 2 term (which in the case of ( |T3| ) is g 2 /X) is crucial 
to the resonance structure of the equation. Since we are 
interested in super- Hubble modes, we need to know for 
which values of the coefficient the mode K = lies in the 
first instability band of the equation. This is the case for 
values 



1 < ^- < 3. 
A 



(16) 



This implies that resonance occurs only in the equation 
( |l3| ) for Q x , and then only for values of g 2 /X in the above 
range (or in the range corresponding to other instabil- 
ity bands). However, if the condition ( |l6| ) is satisfied, 
and if the x field is indeed excited, then parametric am- 
plification of Q x is expected, and via the terms on the 



right hand side of (12), induced exponential growth of 
Qcf, should occur. 

Note, however, that it is not sufficient to show that an 
amplification of Q<p or Q x occurs in order to demonstrate 
that parametric resonance during reheating will have a 
crucial effect on the amplitude of gravitational fluctua- 
tions. Note that also in the conventional treatment of 
fluctuations [EQ] the amplitude of the Q variable grows 
during the interval in which the equation of state of the 
background changes. As emphasized in ||], a straight 
forward way to check if the effect discussed here is a new 
effect is to consider the time evolution of the "traditional 
conserved quantity" £ ( |l9| ) which gives a measure of 
the adiabatic component of the metric fluctuations. In 
the multi-field case, £ is given by |12| 



c 



H 



zQi 



(17) 



In the two field case, the evolution equation for ( is B4] 



*Note that this condition and the corresponding conditions 
in all the other examples discussed in this paper are stable 
against perturbative coupling constant renormalizations, the 
reason being that we consider g 2 ~ A and A should be con- 
strained by the CMB anisotropy results to be a small number 
<< 1. Thus, the perturbative correction terms which are of 
order g 4 are much smaller than either g 2 or A. 



H a 1 



(18) 
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~2 



x 




In the case of a single scalar field, the second term van- 
ishes, but the equation cannot be applied during reheat- 
ing since H = at times when = 0. 

In the single field case, the evolution equation for 
£ which applies also during reheating is (for long- 
wavelength perturbations for which the spatial gradient 
term can be neglected) H 



(l + w)C =0. 



(19) 



This shows that £ is conserved unless w = — 1. Unless 
matter is given by an oscillating scalar field (in which case 
W = — 1 will occur at the turnaround points when (j) — 
0), Eq. (|l^ ) implies that the variable C is conserved on 
scales outside the Hubble radius [jl9| . However, reheating 
corresponds to an oscillating inflaton field, in which case 
the conclusion that Q is constant may break down, as 
discussed in ||. Nevertheless, in the specific single field 
models which have been analyzed jlQ|,|i^p|-|8|] it was found 
that no net increase of Q occurs during the initial stages 
of reheating, and a general proof of the constancy of £ in 
single field models was suggested in || . 

If matter is described, more realistically, in terms of 
multiple scalar fields (each of which is given by a conven- 
tional action), then it appears extremely unlikely that 
w = — 1 will occur at all, since at the points in time 
when = 0, the other fields will not all also be at rest, 
and thus the net value of w will be greater than — 1. 
Therefore, the only realistic possibility for growth of ( is 
as a consequence of the second term on the right hand 
side of ( |I8| ) , a term which corresponds to an isocurvature 
perturbation. Inspection of (|l8|), however, immediately 
shows that during reheating there is the chance of hav- 
ing a very large increase in ( as a consequence of the 
zeros in <fi which arise periodically in time. This effect 
is missed if the scalar fields are treated in the slow-roll 
approximation, or if the change in the equation of state 
during reheating is modelled as a monotonic change from 
a nearly de-Sitter equation of state during inflation to a 
radiative equation of state after reheating. 

However, to see if there is indeed an exponential growth 
of isocurvature perturbations, it is important to take a 
closer look at the equations. In the specific two field 
model of (j^), symmetric initial conditions for the ho- 
mogeneous part of \ giv e X = and x = 0- I n this 
case, it follows from ( p7| ) that £ only depends on Q^,, 
and, since by ( |l2| ) there is no parametric amplification 
of super-Hubble modes of Q</> given that the coupling to 
Q x vanishes, that there will therefore be no paramet- 
ric amplification of Q. The same result can also be seen 
from (|18| ) since for symmetric initial conditions, the time 
derivative on the right hand side of the equation acts 
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FIG. 1. Evolution of (1 + w)( for the model of Eq. (3) 
as a function of the a-dimensional time z = s/\M v \t for 
Xo = Xo = and <f>o = 3.5Af p i , <po = — .1M V \ as initial con- 
ditions for the background. The fluctuation and Q x start 
in the adiabatic vacuum 40 e-foldings before inflation ends. 
The wavenumber is k — 10 2 , which corresponds to five times 
the Hubble radius at the beginning of the simulation. Note 
that the mode is far outside the Hubble radius at the end of 
inflation. 



on a constant. The result is confirmed by our numerical 
analysis (see Figure 1). 

However, due to quantum fluctuations we expect that 
the average of x over a volume corresponding to a partic- 
ular super-Hubble (but sub-horizon) mode will not van- 
ish. It is reasonable [ p5|]2(| to use for the homogeneous 
value of x the r.m.s. value of the renormalized quantum 
fluctuations. It follows from (13) that Q x will experience 
parametric amplification during the initial stages of re- 
heating. It will grow as exp[fMorj\, where fj,o is the Floquct 
index of k = (which for continuity reasons cannot be 
too different from the Floquet exponent for long wave- 
lengths). Via the non- vanishing source terms in (12), 
this will induce parametric growth of Q^, and this quan- 
tity mainly will contribute to the parametric growth of 
Q. Our numerical analysis confirms the above consider- 
ations. In Figure 2 we depict the growth of ( during 
the initial stages of reheating for two different values of 
the homogeneous component of which shows how the 
onset of the parametric growth of £ is dependent on the 
value of x- hi Figure 3 we show how the growth of £ is 
similar to the growth of Q^. 
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z = V\M pl t 

FIG. 2. Evolution in logarithmic 

scale of (1 + w)( for xo = 2 x 10~ 2 A/ P i , xo = (top) and 
Xo = 2 x 10 _8 M p i , xo ~ AMpixo (bottom) as initial value 
for x- The initial condition for cj> and (j> are the same as in 
Fig. 1 in both of the panels. The fluctuations and Q x 
start in the adiabatic vacuum 40 e-foldings before inflation 
ends. The wavenumber is k — 10 2 , which corresponds to five 
times the Hubble radius at the beginning of the simulation. 
The growth of £ is delayed in the second case because the 
background field, and consequently the mixing terms in Eq. 
([t2|), are smaller than in the first case: in this way Q x takes 
longer to feed the growth of Q^, and (. The initial conditions 
for the second case correspond to the values obtained through 
renormalization arguments. 



III. ISOCURVATURE PERTURBATIONS 



Having determined in the previous section that in the 
model (||) there is indeed amplification of £ during re- 
heating, we must now show that this effect is indeed a 
consequence of parametric resonance, and not just an ef- 
fect due to the change in the equation of state between 
the inflationary era and the post-inflationary era, an ef- 
fect which is already taken care of in the "usual" theory 
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FIG. 3. Evolution in logarithmic scale of Q0 (above) and 
Q x (below) for the second set of initial conditions of Fig. 2. 

of isocurvature perturbations in inflationary cosmology, 
which wc dchnc as the results obtained when the tran- 
sition in the equation of state between the inflationary 
slow-rolling phase (p ~ —p) and the post-inflationary 
radiation-dominated phase (p — ^p) is taken to be mono- 
tonic. We will show that in the "usual" analysis there 
can be no exponential increase in the isocurvature per- 
turbation, and that therefore the exponential increase we 
obtain here is a result of parametric resonance. 

The fact that isocurvature perturbations can induce 
an adiabatic component on supcr-Hubble scales has been 
known for a long time |2^- |3l| , |20|| . Entropy perturbations 
act as a source for £ even on scales larger than the Hubble 
radius. This is true even in the case when matter is 
given by a single scalar field. In this case, the evolution 
equation for £ becomes [p2| 



v p dp' 



(20) 



These perturbations, however, are suppressed on scales 
larger than the Hubble radius. 

In models with two or more scalar fields, the equation 
for C is given by (|l^), and it is thus clear that even on 
super-Hubble scales one should expect ( 7^ 0. In the 
approximation in which both fields are slowly rolling, 
the time evolution of £ on scales larger than the Hub- 
ble radius was studied in detail in |24|,|32|,|33| , with partic- 
ular emphasis on calculating the deviations from scale- 



invariance of the resulting power spectrum of density fluc- 
tuations. However, since the analyses made use of the 
slow-rolling approximation, no effects of the dynamics of 
reheating were considered. 

More recently, Taruya and Nambu Jl^] and Bassett et 
al. considered the effect of reheating on the spectrum 
of density fluctuations and discovered a large growth of ( 
due to the initial isocurvature perturbations, however in a 
model in which the necessary initial \ field fluctuations 
are exponentially suppressed during inflation. Bassett 
and Viniegra Jlfj[ then pointed out that the suppression 
would be absent in the model (||). 

It has been known for a long time that isocurvature 
perturbations can be produced in inflationary models 
with more than one scalar field. This issue was initi ally 
considered in the context of axion perturbations in [ p4[ , 
extended to more general two field models in |35| and 
studied in detail in p6fl . It was discovered that initial 
super-Hubble-scale isocurvature perturbations induce an 
adiabatic component by the time that the scales re-enter 
the Hubble radius. 

The gauge invariant expression for the total isocurva- 
ture perturbation in a multi-fluid system is |3l|] 



2 M), 



(21) 



where Spf and Spf are the gauge invariant pressure and 
density perturbations with respect to the total matter 
rest frame and the total speed of sound c 2 s is defined as 
the weighted sum of the i-th speed of sound pi) : 



1 



with 



p E 



= 1 



v., 



(22) 



(23) 



The total isocurvature perturbation can be written as 
the sum of the non-adiabatic pressure component of the 
single component and of the relative isocurvature pertur- 



bation Sij as: 



pr = J2(Sp l - c 2 si 5pi) 



E- 



.(,2 



where 



Si 



_ dpf Sp 



(24) 
(25) 



The relative isocurvature perturbation S%j with respect 
to the total matter frame can be written for our two field 
model as 
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in |36|] in an approximation in which the oscillations of 
the inflaton field were neglected. In this case there is 
no growth of S<p x . An initial isocurvature perturbation 
does induce the growth of an adiabatic component on 
super-Hubble scales, but the final amplitude of the adia- 
batic mode is not much larger than the initial amplitude 
of the isocurvature perturbation, in agreement with the 
earlier analysis in |3J]. 



IV. THREE GOLDEN RULES 
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FIG. 4. Evolution in logarithmic scale of the total 
non-adiabatic pressure pT for the first set of initial conditions 
of Fig. 2. 
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where v is the total perturbed velocity for matter, 
is the homogeneous energy transfer to the <f> component 
{Q<f> + Q x = because of total energy conservation), and 
now Spi are the density perturbations in the longitudinal 
gauge: 



m + V,i 



(27) 



From Equation (p6| ) it follows that the parametric res- 
onance from the matter sector of the theory (the Q vari- 
ables to be specific) induces exponential growth of the 
relative isocurvature perturbation, and hence also of the 
total isocurvature perturbation (see Figure 4). 

In turn, isocurvature perturbations determine the 
change in ( via Eq. ([II]). This shows that in the presence 
of scalar field interaction terms, there is a correlated ex- 
ponential growth of £ and of the relative isocurvature per- 
turbation Srj, x . This exponential growth is a consequence 
of parametric resonance and is absent if the phase tran- 
sition is modelled with a monotonically increasing value 
of w. 

Analogously to Eq. ([l7]) for the Bardeen parameter, 
the total non adiabatic pressure pT can be expressed in 
terms of the Sasaki-Mukhanov variables in the following 
way: 



v ■ v 

—<PiQi - 2V,i Qi - 2 — 

P 3HY, 



2 



(28) 



As mentioned above, the exponential growth of S$ x 
and £ during reheating is a new effect due entirely to 
parametric resonance. The growth of fluctuations in in- 
flationary models with two uncoupled fields were studied 



Based on the analysis of Section 2, it appears that 
several conditions are required in order to have efficient 
parametric resonance of super-Hubble-scale metric fluc- 
tuations. 

1. In the absence of gravitational perturbations there 
must be broad-band parametric resonance in the matter 
sector of the theory corresponding to isocurvature fluc- 
tuations, and k = must be part of the resonance band. 

2. The fluctuations in the matter field which under- 
goes parametric resonance must be effectively massless 
during inflation. More precisely, there should be no large 
net suppression of these fluctuations before the phase of 
parametric resonance. 

3. The homogeneous value of the matter field which 
undergoes resonance must be non- vanishing. This is the 
weakest of the three conditions since it is only required 
if we work strictly to first order in perturbations. 

We show how these three rules are satisfied in another 
model with massless fields, but now based on negative 
coupling instability |p7ft. The potential is the following: 



V{<l>,x) 



\ J.4 2 ±1 2 



l 

T 

r = AA, 



X 



(29) 



with the parameter r = \\ x / g > 1 in order to have a 
potential bounded from below |3^]. This model has an 
attractor for x m the point which minimize the potential 
for X [0 



(30) 



In this way even the third and weakest of the above con- 
ditions is satisfied. The second one is easily satisfied be- 
cause of the negative effective mass for the x fluctuations 
when the background x is small. In order to find the 
unstable bands one can use the attractor solution and 
estimate the frequency of the inflaton and of the fluc- 
tuations Sx during the period of coherent oscillations: 
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If x is small compared to the inflaton </>, then an unstable 
band for k = should be located at 2g 2 = 2A. This gives 
a second order equation for g 2 whose positive root is: 
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FIG. 5. Evolution in logarithmic scale of (1 + w)£ for the 
model @ for 

g 2 = A, A x = 10 2 A (top) and A = X x ,g 2 = A/2(v / 5 - 1) 
(bottom). The initial conditions for the background are 
<t>a = 3.5M p i,0o = - lVXM p i,Xo = 210 -8 x M p i and 
Xo = 0. The fluctuations and Q x start in the adiabatic 
vacuum 40 e-foldings before inflation ends. The wavenum- 
ber is k = 10 2 , which roughly corresponds to five times the 
Hubble radius at the beginning of the simulation. 



2 -1 + y/l + 4A/A x 
9 - A x 7, 



(31) 



We confirm numerically this analytical estimate in Fig- 
ure 5 for two allowed values of g 2 : g 2 = A with A -C A x 
and g 2 = X/2(y/E - 1) with A = A x . The reason for the 
growth of C is similar to the previous case: Q x is para- 
metrically amplified, feeds the growth of Q^, and in this 
case both contribute to the growth of £. 



rolling towards <j> = is terminated by a phase transition 
in the second scalar field x, a field with the double- well 
potential. This implies that during the oscillations of <fi, 
the background value of x is non- vanishing, leading to an 
obvious realization of condition (3) above. 

Parametric resonance in the matter sector of hybrid 
inflation models was studied in detail by Garcia-Bellido 
and Linde f40| , and, in supersymmetric hybrid inflation, 
by Bastero-Gil et al. p] . The resonance of the fluctua- 
tions of the two fields <fi and \ is inefficient for a large set 
of the parameter space since both the 8<j> and 8\ fields 
are effectively massive during the regime of coherent os- 
cillations (Sx becomes massive through the Higgs mech- 
anism). Quite generically, parametric resonance could 
be much more efficient if there is a third field tp which 
couples to both <j> and X- Another interesting possibil- 
ity is to consider a doublet for the field x- Then, even 
in the phase of coherent oscillations, there is a massless 
degree of freedom, namely the "Goldstone" mode. Such 
a situation arises naturally in supergravity models (||] . 

Therefore, as a toy model we will consider the following 
potential for the inflaton field and the doublet x — 



v-2 



v = K^-\x\ 2 ) 2 + \m 2 4> 2 + g 2 4> 2 \xt 

1 2j,2 i 1 2i2( 2 2\ , lvl 

= 2 m(j} 2 9 ^ X2 ' IX 

1 ,/2a 2 i 2\ , ^ 2 2 

- 2 M (X1+X2) + 2X1X2 



xi) 



(32) 



For supersymmetric hybrid inflation, there is only one 
independent coupling constant since 



X= 9 - 



(33) 



The values of the masses and coupling constants are 
constrained by the amplitude of density fluctuations at 
Hubble radius crossing, which is given by 63] 



Sp 
P 
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A 3 / 2 



(34) 



V. A MODEL MOTIVATED BY HYBRID 
INFLATION 

Another natural scenario in which the above condi- 
tions can all be satisfied is hybrid inflation J38|. Hybrid 
inflation is also an attractive framework for implement- 
ing inflation in the context of supergravity models [j39| . 
Since (at least) two fields are involved in the dynamics of 
hybrid inflation, the generation of isocurvature pertur- 
bations is rather natural. In hybrid models, the phase 
of inflation during which the inflaton field <p is slowly 



which should be about 10 5 . In our numerical simula- 
tions, we choose g 2 — 10~ 3 , A = g 2 /2 (as in a supersym- 
metric model), M 2 /M 2 l = 1CT 12 and M 2 /m 2 = 10 10 . 
With these values, the Hubble parameter during the late 
stages of inflation is much larger than m which ensures 
slow rolling of <j). 

In this model, slow-roll inflation takes place while the 
value of the inflaton <f> is larger than <j> c = M / g. For these 
values, the effective square mass m 2 , of x (evaluated at 
X = 0) is positive. Once <f> drops below </> c , the m 2 turns 
negative and inflation ends (at time t c ) via a symmetry 
breaking transition in the matter sector. We choose the 
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basis of the x fields such that the order parameter of 
the transition is \i- However, since in general the initial 
ratio of x fields and x velocities at t c is not the same 



X2 

Xi 



(*c 



Xi 



(35) 



we can with no loss of generality assume that the value 
of X2 at the time t r , when the x transition is complete 
and the oscillations start, does not vanish (as in the 
previous case a reasonable value for xi is the r.m..s. of 
the renormalized quantum fluctuations). Thus, we have 
argued that the third of the conditions mentioned at the 
beginning of this section (non- vanishing background mat- 
ter fields) is naturally satisfied in this model (in contrast 
to the model considered in Section 2). The initial val- 
ues of the matter fields at the beginning of the period of 
oscillation will be 



Xi (t r ) - M/y/X 

0^X 2 (*r) <M/VX. 



(36) 



In contrast to the model (|f|) in which the matter fluc- 
tuations are exponentially suppressed during inflation, 
the effective negative coupling instability in the matter 
fields in the time interval between t c and t r leads to the 
conclusion that the x fluctuations are not suppressed. 
In fact, they are suppressed during slow-roll inflation, 
but build up again exponentially fast in the time interval 
t c < t < t r . To see this, we focus on the evolution of 
the fluctuations in of the x field since, as we shall see be- 
low, these are essential for the effectiveness of parametric 
resonance. We will consider field fluctuations neglecting 
metric perturbations and the mixing terms deriving from 
particle interactions (even if in the hybrid models these 
mixing terms are not perturbatively small and their im- 
portance has been emphasized in plj). Under these ap- 
proximations the evolution equation for $X2 is 



k 2 

5x2 + 3H5x2 = -(^2+9 



Ax 2 



M 



1 5x2 
(37) 



For 4> > <p c , the effective square mass is larger than H 2 
and positive, thus leading to damped oscillatory solutions 



5x2 ~ a 3/2 (t)exp(iujt) 



(38) 



with w = g<f> (in the adiabatic limit). During this time 
interval, however, the homogeneous components of the 
matter fields are also damped. The evolution equation 
for the order parameter xi is 



Xi + 3Hxi 



{g 2 4> 2 + ^xl + ^xl-M 2 ) X i (39) 



Since (up to the contributions from x\ which are negligi- 
ble during inflation) the effective masses in ([39]) and j37| ) 
are the same, the damping rates of xi an d 5x2 are also 
the same for t < t r . In the time interval between t r and 



t r , the signs of the effective square masses in both Equa- 
tion ([39]) for xi an d Equation ( |37| ) for 5x2 are reversed. 
In both cases, the effective m 2 is now —M 2 , leading to ex- 
ponential increase in both xi and 6x2- This period ends 
when xi reaches the minimum of the potential at time 
t r . To summarize the above discussion, the evolution of 
Xi follows 



M 

7a 



Xl(tr 



a M{t r -t c ) 



XI (tc) 



(40) 



o M(t r -t c )-%H(t c -t t ) 



Xi(ti 



where f , is the time at the beginning of inflation and H is 
the Hubble constant during inflation, assumed to be con- 
stant to make the equation simple (this assumption does 
not affect the basic point we are making). In comparison, 
the evolution of 5x2 obeys 



M(t r -t c ) 



5X2 (tc) 



(41) 



z M(t r -t e ) e -%H(t e -U) 



) S X 2(U 



which shows that the exponential growth of 5x2 for 
t c < t < t r precisely makes up for the exponential de- 
cay during the period t\ < t < i c ,exactly as it does 
for the evolution of Xi- Equations ( f40|) and ( flff ) can be 
combined to give 



SX2(tr) 



M/Va 

Xi(U) 



$X2(ti) ■ 



(42) 



This demonstrates that there is no overall suppression 
of the fluctuations in 5x2 before the onset of paramet- 
ric resonance, showing that the second condition for the 
effectiveness of parametric amplification of super-Hubble 
gravitational modes mentioned at the beginning of this 
section is satisfied. 

The final conditions to discuss are the criteria for para- 
metric resonance of the k = modes of the matter per- 
turbations. To do this, we consider the mode equation for 
5x2 during the period in which the inflaton <fi oscillates. 
For general hybrid models, the background dynamics is 
chaotic since both <p and xi oscillate with different fre- 
quencies. However, in the supersymmetric case f4lf| , the 
frequencies coincide and the background dynamics be- 
comes non-chaotic. Both <f> and xi oscillate with the fre- 
quency \[2M. To simplify the analysis, we shall neglect 
the back-reaction of particle production and expansion 
on the inflaton, and neglect the Hubble damping term 
in the equation of motion (this is a good approximation 
since we are considering a case in which H « V2M and 
the fields oscillation are not damped by the expansion of 
the universe). Therefore, we take the inflaton to oscillate 
with amplitude 4> a < <p c , and xi will oscillate about its 
ground state as 



Xi{z) 



M 

71 



(! + /(*))• 



(43) 
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where f(z) is periodic with period 2ir. It is convenient to 
introduce the dimensionless time z = \p2Mt. Denoting 
the derivative with respect to z by a prime, the equation 
for the Fourier mode X2k of 8\2 becomes 



X2k 



X2k( 



,2J,2 



2a 2 M 2 AM 2 



0. 



(44) 



where we have neglected the terms in \2 ■ In the absence 
of the final term (the term containing f(z)), this has the 
form of the Mathieu equation 



X 2 k + X2k (A(k) - 2qcos(2z)) = . 



(45) 



The value of q is q < 1/8, the maximal value being taken 
on if (p a = <f) c , and for long wavelengths A(k) ~ 2q. As 
can be seen from the Floquet instability charts (see e.g. 
Fig. 1 in p7| ) , these values do not correspond to efficient 
resonance. From the evolution of the background fields 
obtained from the full numerical solution of the back- 
ground field equations (see Figure 6) it follows that the 
amplitude of oscillation 4> a is in fact substantially smaller 
than cj) c . In contrast, xi oscillates with a large amplitude. 
Hence, the term containing f(z) in Equation ([li]) is more 
important. This term leads to negative coupling instabil- 
ity (see |37[] for a general discussion of resonant particle 
production by negative coupling instability) for small val- 
ues of k. Hence, we expect parametric amplification of 
long wavelength gravitational fluctuations in our model. 

The above considerations are supported by our numeri- 
cal results. In Figure 6 we show the evolution of the back- 
ground fields (f>, Xii X2 and H as a function of time in 
a simulation with parameters mentioned after Equation 
(|34|), and with initial conditions 4>o — 3M p i, xi — .01Af p i, 
X 2 = .0001M pl . 

As is evident, following an initial transient period the 
three scalar fields oscillate with the same frequency. The 
results for the fluctuation variables Q r f )1 Q X n Qx? an< ^ C 
are shown in Figure 7. The initial perturbation ampli- 
tudes were chosen to be Q<p{to) = 1, Q Xl (to) = Q Xl (to) = ' 
10 -4 , and all their derivatives set to zero for a wavelength 
outside the Hubble radius (k = 0). As is evident, after 0, 
the initial transient period, all four quantities grow al- 
most with the same Floquet exponent, as expected from 
our analytical analysis. 

At this point, an obvious question is whether the field - 
X2 is essential in order to obtain parametric resonance 
of super-Hubble-scale cosmological fluctuations. In fact, 
the field fluctuations 5xi also will experience an effective 
negative coupling instability , and therefore the pres- 
ence of X2 is not essential for this supersymmetric choice 
of the parameters. The equation of motion for 5xi for a 
hybrid model with two field is 

,k 2 

,2 




H 



0001935 



0001925 
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Sxi + 3H5xi 



(-+ 5 V + 3A X ?-M 2 )<5 X1 . (46) 



z = Alt 

Evolution of the background dynamics for a su- 
model of hybrid inflation with g 2 — 10 -3 , 
= 10" 12 and M 2 /m 2 = 10 10 . The initial condi- 
■o = 3M p i, xi = -OlMpi, xa = .OOOlMpi and all the 



FIG 
persymmetric 
M 2 /M 2 vl 
tions are 
field velocities set to zero. 



The effective squared mass is large and positive during 
slow-rolling. At the beginning of the transient period 
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FIG. 7. Evolution in logarithmic scale for the fluctua- 
tion variables Q^, Q X1 , Q X2 and (1 + w)(. The initial 
perturbation amplitudes were chosen to be Q<f,(to) = 1, 
Qxi (to) = Qxi (to) = 10 _4 M p i, and all their derivatives set to 
zero for a wavelength outside the Hubble radius (k — 0). 



(when xi starts rolling down its potential but is not 
yet close to the minimum of the potential) the effec- 
tive squared mass turns negative. Once xi gets close to 
its equilibrium position, the effective squared mass will 
again be large and positive (the factor 3 in the third term 
on the r.h.s. of (|4q ) is crucial). However, since Xi is os- 
cillating with a large amplitude, the effect of the large 
mass will be periodically cancelled out by these oscilla- 
tions. Neglecting the expansion of the background, Eq. 
(Eq) can be written as 



Xlk 



Xik( 



2a2 



2a 2 M 2 2M 2 



0. 

(47) 



From our numerical results (Figure 6) we expect that 
the amplitude of f(z) will be only slightly smaller than 
1. Hence, we expect negative coupling instability for long 
wavelength metric perturbations also in the two field case 
|fli|| . However, as demonstrated below for values of the 
coupling constants which do not correspond to the su- 
persymmetric point, the Floquet index in the two field 
case will be smaller than in the three field model. 

Other interesting effects happen if we go away from the 
supersymmetric point g 2 — 2A. In analogy with the re- 
sults of f43| which show that random noise in the inflaton 
leads to an increase in the strength of the parametric in- 
stability, we expect that the chaotic background dynam- 
ics will not eliminate but rather strengthen the resonance. 
Chaotic background dynamics is expected for g 2 ~ A in 
hybrid models [ ^0[ . This effect is shown by our numerical 
simulations (Figure 8 and Figure 9) which show that the 
parametric resonance of super- Hubble-scale gravitational 
fluctuations for the choice g 2 = A is larger than in the 
supersymmetric case, where no chaoticity is present [ fi"l|| . 
Figure 10 shows how the presence of the "Goldstonc" 
mode xi changes the development of the resonance in 
this chaotic case. In the two field case the Floquet in- 
dex with which £ grows is smaller than the corresponding 
index in the three field case. 



VI. DISCUSSION 

We have studied the parametric amplification of long 
wavelength gravitational fluctuations during reheating in 
two field inflationary Universe models. We have par- 
tially confirmed the results of Bassett and Viniegra |l(t | 
and shown that this effect is possible for certain models. 
We have established criteria under which an exponential 
increase in the amplitude of cosmological perturbations 
during the period when the inflaton field oscillates should 
be expected. It is crucial that there must be either broad- 
band parametric instability or negative coupling instabil- 
ity in the matter sector of the theory (i.e., in the absence 
of gravitational perturbations). This will excite isocurva- 
ture fluctuations during reheating. It is important that 
there be no net exponential damping in the amplitude of 
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FIG. 8. Evolution of the background dynamics for a model FIG - 9 - Evolution in logarithmic scale for the fluctuation 

of hybrid inflation with g 2 = A. The parameters and the variables Q , Q X1 , Q X2 and (1 + w){ for the case g = A of 
initial conditions for the fields are the same of Fig. 6. Fi S- 8 - The initial perturbation amplitudes are the same of 

Fig. 7. 
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z = Mt 



FIG. 10. Evolution in logarithmic scale of (1 + w)( for the 
choice g 2 — X for a two field hybrid model (x2 = 8x2 = 0). 
The parameters, initial conditions for the background and for 
the perturbation amplitudes are the same of Figs. 8 and 9. 



the isocurvature fluctuations before reheating. The reso- 
nance in the matter sector then induces a resonance in the 
gravitational sector provided that the background values 
of the matter fields do not vanish. Since large coupling 
constants are not necessary in order to have efficient res- 
onance, the effect is stable against perturbative coupling 
constant renormalizations. Wc have shown that in this 
case the resulting increase in the amplitude of the adia- 
batic mode, conveniently tracked in terms of the variable 
C, and of isocurvature fluctuations, tracked in terms of 
the non adiabatic pressure pT, is exponential and is due 
to the oscillations in the inflaton field. This means that 
the effect is absent if the phase transition is modelled by 
a monotonic change in w = pj p. 

We then argue that the conditions under which para- 
metric amplification of long wavelength gravitational 
fluctuations occurs are naturally satisfied in a class of 
models of hybrid inflation. The presence of a complex 
matter scalar field enhances the resonance, since it en- 
sures the existence of a field which is massless in the 
true vacuum of the theory, but it is not crucial if there 
is negative coupling instability. However, note that the 
existence of massless modes is helpful for the effect to 
occur. Such massless modes arise quite generically in 
string theory (see e.g. for a recent review). Thus, 
the parametric amplification of long wavelength fluctua- 
tions may be also present in models of inflation based on 
string theory. 
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